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Abstract 

We calculate inter-meson potentials numerically by solving classical equations 
of motion derived from Dual Ginzburg-Landau (DGL) Theory. Inter-meson 
potentials in DGL theory are shown to be similar to those of the string- 
flip model and well reproduce behaviors of the short-range interaction at the 
classical level. We also compare our results with those from lattice QCD 
Monte carlo calculations. 
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I. INTRODUCTION 



Quark confinement is one of tlie main ploblems tliat sliould be solved in liadronic pliysics. 
Tlie linear potential between quark and anti-quark at long distance can be calculated by 
Monte Carlo calculations of lattice QCD and the non-relativistic quark potential model 
explain the low-lying hadron spectrum with the infinitely rising potentials. [jT| 

Quark potentials in multi-hadron system should also be investigated. A difficulty appears 
when one try to understand multi-hadron systems in the framework of linear rising quark 
potential model. It gives a long-range attraction, called color van der Waals force, between 
two color-singlet hadrons which seems to contradict experimental data of nucleon-nucleon 
scattering. |^ 

In treating multi-hadron system, the string flip model was proposed phenomenologically. 
In the string flip model, hadrons do not interact each other except when the strings flip 
into another combination. Hence the string flip model is by construction free from color van 
der Waals force. But it is not so straightforward to understand the string flip model starting 
from QCD. 

Dual Ginzburg-Landau (DGL) theory is an infrared effective theory of quark conflnement 
derived from QCD. |^ 'tHooft and Mandelstam conjectured that QCD vacuum may be dual 
to superconductor and that monopoles play the analogous role of Cooper pairs in QCD. ||5|,|^ 
'tHooft proposed abelian projection to extract relevant dynamical variable at low energy in 
QCD and suggested that the variable is color magnetic monopole. The abelian projection 
is a prescription of partial gauge flxing which reduces SU{3) gauge group to U{1) x U{1) 
and monopoles appear as a point like singularity. Following Bardacki and Samuel we can 
integrate out monopole trajectories. [§] Introducing dual vector gauge flelds, we can 
get Dual Ginzburg-Landau Lagrangian that is an abelian effective theory with magnetic 
monopoles. 

Meson, baryon and spin-dependent baryon potentials have been calculated numerically in 
DGL theory. pO|JTl[| The model also explains the characteristic features of flnite-temperature 
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transition of pure QCD found by lattice calculations. ||T^ Monopole condensation and chiral 



symmetry breaking is also discussed. 13 



In this paper, we treat multi-hadron system (for simplicity meson-meson system) in DGL 
theory. We numerically solve classical equations of motion derived from DGL theory and 
obtain static potentials V2 in two body systems and V4 in four body systems. Inter-meson 
potential is evaluated from E4 = V4 — 2V2. We also compare our results from DGL theory 



with those from lattice Monte Carlo calculations done by Green et al.. [|14 



II. INTER-MESON POTENTIAL IN DGL THEORY 

DGL theory is represented by a Lagrangian 

Cdgl = -j{H^^) +Y\\{d^ + igta- C^,)Xa - >^[\Xa? - v"^) \ (1) 

^ a=l ^ ^ 

with the constraint Yl\=i ai'g(Xa) = , where = (C^, C®) is the dual abelian vector 
field with respect to f/(l)xf/(l) and Xa{ci = 1,2,3) are the complex scalar monopole fields 
which couple to covariantly under magnetic f/(l)x[/(l). Here are the root vectors 
: ei = (1,0), €2 = (-1/2,-^3/2), eg = (-1/2,^3/2). H^, is defined as 

H^.. = d^C, - d,C^ + e^^A. / A n^{n ■ d)-\x - y)feAy). (2) 

where is an arbitrary constant four vector and jgx = Us^js) is an external color-electric 
current. We have already assumed monopole condensation and abelian dominance 0] which 
is supported by Monte Carlo simulations of QCD. ||l5l It is noted that this form is just dual 
to relativistic Ginzburg-Landau theory with an external magnetic current. 

We consider for simplicity the following symmetric configurations of four static quarks 
in three dimensional space : 



-^(x).(.-f).(. + f) 
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+ + + (3) 

where 

4 = (QLg^) = (-|,-^). (4) 

Here color-electric coupling e is related to g by Dirac quantization condition eg = Att. 

We now look for static solutions in which the time components are neglected. Adopting 
a unitary gauge aig{xa) = , we get equations of motion : 

3 

V X E~2g^Y.^o.{ea-C)\xaf = 0, (5) 

a=l 

A\xa\ - 2g\e^ ■ C)\xa\ - 2X{\xa\' - v')\xa\ = 0, (« = 1, 2, 3). (6) 
The energy density £ is written as 

^ = \E+j: \{V\x.\f + g\e^ ■ Cyix^f + X{\Xaf- v^y] . (7) 
Here the color-electric field is defined as 

E = V xC + Es, (8) 

— * 

where Es is the singular string field which vanishes everywhere except on the Dirac strings. 
It comes from the last term of Eq. (Q). 

If we require E ^ and Xa ^ v at the spatial infinity, our gauge choice makes the 
C field vanish at this region as seen from Eq. (|^). This means that the string singularity 
does not extend to infinity but connects the quark charges in the unitary gauge. There are 
various ways to connect the charges. 

Here we consider the Dirac quantization condition which guarantees unobservability of 
the string. We can change the string location by a single-valued gauge transformation, if 
the integral of the dual vector field along an infinitesimal closed loop is quantized as follows 
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\imjCdl = Y.i,Qi , (9) 



where ^ is an integer and eg = An is used|^ Condition (|^) selects two-type configurations of 
tfie Dirac string (Figs. 1 (a) and (b)) out of those which connect the charges (Figs. 1 (a), 
(b) and (c)). Note that in the unitary gauge the location of the color-electric flux coincides 
with that of the Dirac string where the value of the monopole fields Xa vanishes. 

From = y/Sjs it is expected that there exists a solution satisfying = V^C^ = C^. 
Then two equations of (^ (a = 1, 2) become identical, so that we get the same solution for 
Xi and X2 (|xi| = 1X2 1 = x) ■ Moreover we can solve the remaining equation of (||) (a = 3) 
with the condition that the static energy is minimized : 

|X3|=^^. (10) 

To solve Eqs. (|^) and (^) numerically, we have to avoid infinite quantities such as the string 
singularities. Following Ball and Caticha, we separate 



C = Cd + C. (11) 

Here Cd represents the Coulomb field generated by the four charges. In the case of (a) type 
string location in Fig. 1, it is written as 



Ql cos Oi + l 2 Qi COS 612 + 1 



>i — 1 r-. — TTH^^ 



47r I r — ri I sin 9i 47r | r — r2 1 sin 62 

Ql COS^3-|-l Ql COSe^ + l ^ 

An \r — \ sm U-^ An \r — sm 64^ 



where the notation is explained in Fig. 2 . The new fields C and V x G are well behaved 
everywhere. The field equations are reduced to 



^Consider an infinitesimal closed loop around a Dirac string. After a gauge transformation around 
the closed loop, the dual vector gauge fields C = C + V Am become regular : lim §{C + ^Am)dl = 
0. Also the Xa fields remain single-valued : gea'{§ VAm^O = 27rn. Hence we see gea-i).vai§ Cdl) = 
27rn which is the generalized Dirac quantization condition m-q = 27rn. Using Eq. (^), we get Eq. 
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V\ - g\Cn + Cfx = 2X{x' - (13) 
V'C-V- (VC) = 3(?V(Cd + C'). (14) 

Eq. ([TlD not only eliminates the singularities from the equations for C but also enables 
a separation of the infinite Coulomb self-energy. The static potential energy in four body 
system without the self-energy is given by 



+ J Sx {-23\^C ■ C + 2A(i.' - x')} . (15) 

As the boundary condition at the space-like infinity, we have required C — * and % ^ f . 
This means the vacuum of the dual superconductor is realized at the infinite boundary. 

The static potential V2{R) in two body system can be evaluated in a similar way. Finally 
we get inter-meson potential i?4 using the equation 

, 2V2 (i?), \iR<D, 
E^{R,D) = V^{R,D)-{ (16) 

2V2iD), ■iiR>D. 

III. NUMERICAL METHOD 

As we treat the system numerically, we rescale the variables as = fx, C* = 
v~^C, C'^ = v-^Cd, x' = v-^x and V\R, D) = v-^V{R, D) . 

The spatial boundary conditions are imposed on the surface of sufficiently large but finite 
volume. Since the system has (anti) symmetries under x — — x, y —y or z —z , we 
can restrict ourselves to the region (x>0, y > 0, z > ) with the appropriate continuity 
conditions. The region is discretized into an Nx x Ny x Nz lattice. Since a fine lattice is not 
needed near the boundary surface, we have used lattices whose lattice spacing becomes larger 
as the distance from the quark sources. We have used the Gauss-Seidel method applying the 
successive overrelaxation technique. Furthermore, we have adopted the following procedure 
to get rapid convergence. We work first on a coarse lattice and next employ the solution as 
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the intial configuration on a finer lattice which is constructed by putting new sites between 
the original lattice sites. And we repeat this process. The largest values of N^, Ny, and 
attained are 80, 120 and 120, respectively. 

The procedure of changing the lattice size shows us the lattice spacing dependence of 
the energy V4 and V2. It is nicely fitted by ¥4^2 = Cq + Cia^ + C2a^ where a is a typical lattice 
spacing. This dependence gives the continuum limit of the energy. 

Varying the inter-quark distances R and D, we get the potential V4 and V2. 



We first illustrate the results of the inter-quark potential. As shown in Fig. 3, inter-quark 
potential in DGL theory is Yukawa-like in the short distance and linear in the long distance. 
DGL theory realizes quark confinement. 

To compare our results based on DGL theory with those from SU{2) lattice QCD cal- 
culation obtained by Green et al., |T^ we study the SU{2) case. In SU{2) case, Eqs. (|T3|) 
and dU) are almost the same as in SU (3) case except that the coupling 3g^ is replaced by 
2g^ in Eq. ([HI). The four quark potential is modified into 



Our results in SU{3) case is qualitatively similar to those in SU{2) case. To compare these 
results, we have to adjust the parameters in both results. Free parameters are e. A, v in 
DGL theory and 02.4 which is the lattice spacing at /3 = 2.4 in SU{2) lattice QCD. Inter- 
quark potential V2{R) calculated numerically in DGL theory gives the scaled string tension 
a/v^ = 6.2809(2), whereas lattice SU{2) QCD Monte Carlo simulation gives the scaled 
string tension 0"a2,4 = 0.0699. [|T^ The string tension a = (440 MeV)^ determines v and 
02.4 to be 175.567(3) MeV and 0.1179 fm, respectively. Because all charged gluons as well as 
neutral gluons contribute equally in the short distance, the Coulomb term of V2{R) in DGL 



IV. RESULTS AND DISCUSSION 




(17) 
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theory should be equal to one third of the term — 0.249/i? derived from the lattice data 



Hence we set e = 2.04. In Ref. [O], it is found that the QCD vacuum is near the border 



between type 1 and 2. Hence we take A = |(^)^ in accordance with the Ginzburg-Landau 
parameter k = VX/ g = 1/ V^- 

Now we can calculate the inter-meson potential. We consider both cases of Dirac string 
location (Fig. 1 (a) and (b)). Using Eq. (|16|) we evaluate the inter-meson potential Ei{R, D) 
in both cases and choose the lower one. In Fig. 4 the results are shown when R is fixed. Two 
mesons hardly interact each other in the region D ^ R. But as two mesons are getting closer, 
interaction between them becomes sizable at some value of D. The maximally interacting 
point is just the square {R = D) point. Still decreasing D makes Dirac strings flip to take 
the holizontal direction and the interaction gets weaker. This behaviour is similar to the 
ansatz made in the string flip model. 

To see whether these inter-meson potentials in DGL theory are long-range or short-range, 
E^[R, D) in the region D > R is fitted by the function 

p-SD 

- a-lT- (18) 

when R = 0.1179 fm. This function reproduces well the data when 6 = 1.43(1) GeV and 
7 = 1.145(1). So the inter-meson potential in DGL theory is short-range. 

When R is taken larger, the peak value of interaction becomes smaller (Fig. 5). The case 
R = D is shown in Fig. 6. In Fig. 7 we also show the potential in the case that Dirac strings 
run in the same direction. In this case the inter-meson potential is repulsive as expected. 

We also compare our results with those from lattice QCD . |jl^ Both results are similar 



qualitatively as seen from Figs. 4, 5 and 6. The degree of falling at i? ~ D in DGL theory 
is a little weaker than that in lattice QCD. This difference may come from the assumption 
that we have neglected the off-diagonal gluons whose contribution is sizable in the short 
distance. 

In Fig. 8 the color electric flux distribution is shown. The color electric flux spreads in 
all direction in Coulomb phase (f = 0), while it is squeezed between quarks which belong to 

8 



the same meson in the confinement phase {v ^ 0). This illustrates the interactions between 
mesons are short-range. Thus DGL theory reproduces both properties of quark confinement 
and asymptotic separability of mesons at the classical level. The potentials in DGL theory 
is linear in the color electric fiux direction, whereas Yukawa-like short-range in the direction 
perpendicular to the fiux. 

Recently an attempt to determine the couplings in the abelian effective monopole action 
has been performed. Monopole action can be derived numerically from the monopole current 
configurations given after abelian projection on the lattice. This action is equivalent 
to non-compact abelian Higgs action on the dual lattice in the long-range region. ||19[ It is 
very interesting to fix the couplings in DGL theory directly from QCD. 



This calculation were performed on Fujitsu VPP500 at the institute of Physical and 
Chemical Research (RIKEN). 
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FIGURES 
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(a) 



(b) 



(c) 



(d) 



FIG. 1. Locations of Dirac strings, (a), (b) each charge has one Dirac string. Two Dirac 
strings run in the opposite direction, (c) each charge has two Dirac strings, (d) each charge has 
one Dirac string. Two Dirac strings run in the same direction. 




FIG. 2. Coordinates used in Eq. (IT 
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FIG. 3. Inter-quark potential derived from SU{2) DGL theory when e=2.04, k = 0.707, 
cr = (440 MeV)^ 
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FIG. 4. Inter-meson potential dreived from SU{2) DGL theory (dot) when e=2.04, k=0.707, 
a = (440 MeV)^. The open circle denotes potential from SU{2) lattice QCD in Ref. |14]. Quark 



charges are located at the corner of the rectangle and two Dirac strings run in the opposite direction 
(Figs. 1 (a) and (b)). R is fixed to be 0.1179 fm. 
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FIG. 5. Inter-meson potential dreived from SU(2) DGL theory when R is fixed to be 0.2358 fm. 
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R = D [fm] 



FIG. 6. Inter-meson potential dreived from SU{2) DGL theory in the R = D case. 
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FIG. 7. Inter-meson potential dreived from SU (2) DGL theory when two Dirac sitings run in 
the same direction (Fig. 1 (d)). R is fixed to be 0.1179 fm. 
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FIG. 8. Color electric flux in confinement phase v 7^ (solid line) and that in Coulomb 
phase V = (dashed line). Color electric charges are located at (x, y)=(±0.795 fm, ± 0.445 fm). 
Parameters are taken to be e=2.04, k=0.707, a = (440 MeV)^ . 
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